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Abstract. 

The quantum ANNNI chain in a transverse field is investigated by means of 
the bosonization approach in the limit of large next-nearest neighbor interaction. 
In this regime, this model can be viewed as a weakly coupled two-leg zigzag ladder 
which enables us to derive its low energy effective field theory. In particular, it is 
shown that the effect of frustration in the system is captured by the presence of 
a non zero conformal spin perturbation that accounts for the existence of all the 
incommensurate phases of the model. 
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One of the most striking effect of frustration in magnetic systems is that it can lead to a 
huge number of degenerate ground states. As a result, the nature of finite temperature 
phases may not, in contrast with ferromagnetic systems, be linked only to the sole 
nature of the microscopic degrees of freedom and the dimension of space-time. Since 
frustration induces strong fluctuations that involve large number of spins, the low 
energy physics is determined by interferences between very different regions of phase 
space. With this picture in mind one may not expect the field theoretical description 
of the frustration to be an easy task. Of course there have been several attempts to 
describe frustrated magnets by field theories. But, to our knowledge they were mostly 
restricted to models that displayed well defined ground states, such as helical ordering, 
and the main effect of frustration was captured by an enlargement of the dimension 
of the order parameter space jjj. Here our aim is to single out an operator, in the 
Renormalization Group (RG) sense, that captures the effect of frustration. In this 
respect, we shall consider the most studied frustrated system i.e. the two dimensional 
ANNNI model (see Refs.[^[ ||, [| for a review). This model is characterized by a 
nearest-neighbor ferromagnetic interaction (— Ji < 0) and a competiting next-nearest 
neighbor antiferromagnetic interaction (J2 > 0): 



where denotes nearest neighbor pairs and [k,l] next-nearest ones along a single 
space direction. The phase diagram of the model (|) is very rich and well known and 
can be summarize as follows. At zero temperature, the ground state is ferromagnetic 
for J 2 < Ji/2 and antiferromagnetic (actually an antiphase TtJ.4) f° r J2 > Ji/2- At 
the special point J 2 = the ground state is infinitely degenerate and correlation 

functions are short ranged. At finite temperatures (T ^ 0), one observes five different 
phases depending on the parameters (see figure (1)): ferromagnetic (F), paramagnetic 
commensurate (PC), paramagnetic incommensurate (PI), incommensurate critical 
phase (IC, also usually called "floating phase") and antiphase (A). Three different 
kind of phase transitions are present in this diagram: An Ising transition between the 
F and PC phases, a commensurate-incommensurate transition^ between the A and 
IC phases, and a Berezinski-Kosterlitz-Thouless (BKT) transition separates the IC and 
PI phases. There is also a disorder line that extends down to zero temperature which 
divides the PC and PI phases^]. As a most notable feature of this phase diagram is 
the existence of an incommensurate critical phase in a finite region of the parameter 
space. In a seminal work, Villain and Bak Q| have proposed an approximate effective 
theory in terms of fermions valid in the vicinity of the degenerate point J2 = Ji/2. All 
their predictions have been further confirmed by numerical investigations and series 
expansions. It is the purpose of the present letter to propose a complementary and 
alternative low energy description of the ANNNI model valid in the large J2 limit. 
Furthermore, within our approach, we exhibit a particular operator that is at the 
origin of the incommensurate phases found in the phase diagram (figure (1)). 

Our starting point is to map ([!]) into a quantum Ising model in a transverse 
magnetic field. This is done by introducing an anisotropy between the nearest neighbor 
interaction in the ^-direction (Ji) and in the y-direction ( J T ). As it is well known (see 
Ref. Q for a review), the physics of the classical model (|l|) should be equivalent to 
the one described by the following quantum Hamiltonian: 
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Figure 1. Phase diagram of the classical two-dimensional ANNNI model. 

where c^, cx^ are Pauli matrices and 2/3 J* = lncoth(/3J T ). The model (||) is nothing 
but the one-dimensional axial next-nearest neighbor Ising (ANNNI) model in a 
transverse field. As seen in figure (^), it can also be viewed as two quantum Ising 
chains labelled (1) and (2) coupled by a "zigzag" interaction with strength J\ and 
thus identifies with a two-leg quantum Ising ladder with Hamiltonian: 

U = E [-!3J> Z a(n+l)-0J2O-*(n+lK(n+l + l) 

n,a— 1,2 

- E CT i (» + \) [ ff 2 (») ~ °f (»» + 1)] + (1 - 2). (3) 

n 

Notice that in order to obtain (^) we have performed an (unphysical) gauge 
transformation on the &th chain (a = 1,2): er§(rc + a/2) — * (— l)" +a cr^(n + a/2). 
As we shall now see, the model (||) can be consistently described by a continuous field 
theory in the limit J± << J2 and J2 ~ J*. 



-J 2 
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Figure 2. In the large J2 limit, the ANNNI model is better seen as two weakly 
coupled quantum Ising chains. 

Continuum limit. We shall study the model in the vicinity of the 
antiferromagnetic phase where J\ << J2- One can take advantage that in the limit 
J% = the model is equivalent to two decoupled Ising models which are critical when 
J2 = J*. One can therefore expand the theory around the conformal invariant fixed 
point with symmetry Z2® Z2. In the critical regime the two Ising chains are described 
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by two pairs of right and left Majorana (real) fermions ift a (R,L)i a — 1,2: 

7~to = ^2 (^aRdxtpaR ~ ^ahdx^ah) ~ TH (e% + £ 2 ), (4) 

0=1,2 

e a = iipaRipaL being the energy operator of the &th Ising model whereas the mass gap 
and the velocity are given by: m = 2(J* — J%)jT « 1, v = a J*/T (ao being the 
lattice spacing). In the absence of interchain interaction, the system is disordered 
(respectively ordered) when m > (respectively m < 0) and one has (a a ) = 
(respectively (a a ) ^ 0) where a a are the order operators associated with the two 
Ising models. We consider now the interacting case and take the continum limit of 
the second term in equation (^) in the regime J\ -C J 2 to obtain: 

H~H +gV, V = al /4 (aid x a 2 - o- 2 d x o-i) . (5) 

When m ~ and g = J\/2T = 0, the model is conformal invariant with central 
charge c = i + h = 1. In the generic case, this fixed point is perturbed by the 
thermal operators e 1 ( 2 ) which have scaling dimension A = 1 and the operator V which 
has scaling dimension A = 5/4. It is important to stress that both operators reflect 
very different physical behaviors. While the former is the standard operator measuring 
deviation from criticality in non-frustrated Ising magnets, the latter encodes the whole 
effect of frustration and is responsible of the non trivial phases of the ANNNI model as 
we shall see. The operator V manifests itself by the fact that it is a parity symmetry 
breaking perturbation with a conformal spin equals to 1. The effect of such non- zero 
conformal spin term is non trivial since the usual irrelevant versus relevant criterion 
does not hold for such a non-scalar contribution (see for instance the discussion in Ref. 
[||). In fact, a similar operator also appears in the study of the S=l/2 Heisenberg 
zigzag ladder Q and has been called a "twist term". Such contribution is difficult to 
handle non-perturbatively but it has been stressed that this term may represent a 
new mechanism for incommensuration in one-dimensional systems^, |To| . As we shall 
see, in the particular case of the ANNNI model, the effect of this operator can be 
elucidated by means of the bosonization approach and accounts for the formation of 
the non trivial incommensurate phases depicted in figure (1). 

Bosonization and Effective Theory. We shall use the well known equivalence 
between two critical Ising models, characterized by four chiral real fermions "0ifl(L) an d 
i ) 2R( l) , an d a free boson theory described by the chiral fields 4>r(l) ■ The bosonization 
rules are (see Ref. || for a review): 

(tpi + iip 2 )R(L) = —X= exp(±iV4n<j)R(L)), (6) 



from which it follows that Tto is equivalent to a sine-Gordon model at 1 = Air (free 
fermion point) . The bosonization of the twist term (|^) requires more work. A suitable 
bosonic expression is obtained by considering the following operator product expansion 
(OPE) that stems from the fact that the order operator cr a is a primary field: 

(Ti - T 2 ) (z)aia 2 (w, w) — — - — \a 2 dax - o\do 2 \ (w, w), (7) 

z — w 

where T a {a = 1,2) are the energy momentum tensors in the holomorphic sector 
(z = vt + ix) associated with the Ising models. The fields in the left hand side of 
equation (^) can be expressed in terms of the bosonic fields: T\ — T 2 = cos(\/167r0i) 
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and oi<72 = \/2 sin(v / 7r<i>) ($ = fa + 4>r). We thus deduce by performing the OPE in 
the bosonic theory the following representation of the twist term: 

n = v - ((a,$) 2 + (c^e) 2 ) - cos - \ 9 -^- cos sin (Vi^e) , (8) 

Q = fa — fai being the dual field. The next step of our approach is to map (||) onto 
the XXZ Heisenberg model in a magnetic field: 

H = ^ (4 + J|) + 4™ (J£J£ + Jp v L ) + irgV2(jy - Jl) + 2ttAo J z R Jh (9) 

where u — 5u/4, Ao = 3v/2 and the SU(2)i Kac-Moody currents are given by 
J£ R = dxfa,rt/V2n and = exp (±iv / 87r</>L i _R)/(27rao). The two expressions (g) 
and (^|) can be shown to be equivalent by a canonical transformation at the special 
value Ao = 3v/2. As a result, in an appropriate basis, the twist operator acts as a 
magnetic field. The model (0) is more conveniently analyzed in a basis where the 
magnetic field lies along the z axis. To do so, we perform a 7r/2 rotation around the 
x-axis in the spin space so that one finally obtains after rebosonizing once again the 
currents: 

H = Y + ( 9 - e ) 2 ) - 9i cos (v/8^Q$) - 92 cos (v/87r/Qe) - h d x O, (10) 

where 



/ 1 — mao /U — : ; : 

Q = \ — 1-, u =u^{l + ma /u){l-ma /u), (11) 

y 1 + ma /u 

9i(2) = (2m ± Ao/ao)/(47rao), and ft, = g^/n/Q. The effective field theory ( |To| ) is the 
main result of this work and all the different phases (apart from the F and PC phases) 
observed in the ANNNI model can be deduced from a simple analysis of it. 

Phase Diagram. Consider first the high-temperature phase, i.e. when m > 0. 
Since Q < 1, for sufficiently large to, the term cos(y^7r/Q9) is strongly irrelevant 
and can be dropped. The remaining theory is then easy to analyse and a mass gap 
to all excitations is generated due to the presence of the relevant cos(v / 87rQ$) term. 
On the other hand the operator d x Q leads to incommensurate fluctuations of the 
O dependent correlation functions. This phase corresponds to the incommensurate 
paramagnetic phase PI. This picture is confirmed by an exact solution at the special 
value Q — 1/2 where the model becomes equivalent to that of free fermions with 

dispersion e\.(k) = v 2 (k± g\f2ir/v) + 9w 2 /16aQ. As readily seen, there is a spectral 
gap and incommensuration develops as soon as g ^ with wavevector ft* = y/2ng/v. 

Similarily one can study the low temperature regime i.e. m < and Q > 1. 
For sufficiently large to, it is now the cos(v / 87rQ"I>) term that is irrelevant and can 
be dropped. After a duality transformation, the resulting Hamiltonian is equivalent 
to the XXZ chain in a magnetic field along the z axis which is equivalent to the 
bosonized version of the Villain-Bak theory ||TTf derived in the vicinity of the degenerate 
point J2/J1 = 1/2, T = 0. At small g, there is a gap to all excitations with 
no incommensuration: it corresponds to the A phase. As g grows, the magnetic 
field increases until it reaches the gap at some critical value (g c i) above which 
the excitations become massless. For g > <? cl , the system displays as well an 
incommensurate behavior with wavevector k* ~ \/ g 1 — g 2 c \ and one enters the floating 
phase IC. The nature of the transition at g — g c \ is of a commensurate- incommensurate 
tyP e Pl- Finally, as g further increases, the cos(v / 87rQ$) term eventually becomes 
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relevant and opens a gap at a critical value g C 2 where a BKT transition to the PI 
phase takes place. This picture is, as above, confirmed by an exact solution at the 
point Q — 2 where the model ( [To| ) becomes equivalent to massive free fermions. 

As seen, the previous analysis correctly reproduces the phase diagram of the 
ANNNI model in the vicinity of the A phase where frustration plays its tricks. To 
do so, we have assumed that \m\ was sufficiently large to be able to neglect one of 
the cosine terms in equation (|l^) . One can question the validity of this scheme when 
m ~ where both cos(v / 87rQ"I>) and cos(-\/87r / Q&) operators, being almost marginal, 
compete. A detailed analysis of the RG equations associated with (|l(j) is thus called for 
but one is faced with the difficulty that the coupling constants are not small since Ao 
is of order one. One has therefore to make the hypothesis that the qualitative feature 
of the RG approach does not depend on the strength of g\ and and treat them as 
small couplings. The RG equations have already been obtained in a different context 
by Giamarchi and Schulz|l2]|. It follows from jl2| that the results we have obtained for 
large \m\ remain valid for small m confirming our hypothesis. However, there is room 
for a spin-flop transition from the A phase directly to the PI phase if the cos(^/87rQ$) 
perturbation blows up before the d x Q term closes the gap. The occurence of such 
a transition strongly depends on the microscopic couplings of the bare Hamiltonian. 
The analysis of the RG equations for Ao = 6m/5 reveals that a spin flop transition does 
not occur for the ANNNI model. However, care has to be payed since for this value 
of Ao , perturbation theory strictly does not apply and this leaves open the question of 
the existence of a Lifshitz point in the ANNNI model. Furthermore, notice that the 
presence of a Lifshitz point is a generic feature of ( |l0| ) . Therefore one may expect that 
other lattice Hamiltonians displaying the same qualitative phase diagram than the 
ANNNI model will exhibit such a point. In summary, we have derived a low energy 
description of the ANNNI model in the large Ji limit where the system can be viewed 
as a weakly coupled two-leg zigzag ladder. This enables us to start from a conformal 
invariant fixed point in the vicinity of which the continuum limit is well defined and 
frustration manifests itself through the twist operator. This approach accounts for all 
the incommensurate phases of the ANNNI model and in the low temperature limit 
matches the Villain-Bak description near the degenerate point. 

D. Allen is a CNRS post-doctoral fellow and thanks l'universite Cergy-Pontoisc 
for hospitality while this work was initiated. 
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